Critical shell thickness for lnAs-AI ;i lni_ x .As(P) core-shell nanowires 

C. M. Haapamaki, 1 J. Baugh, 2 ' 3 and R. R. LaPierre la) 

^Department of Engineering Physics, Centre for Emerging Device Technologies, 
McMaster University, Hamilton, Ontario, L8S 4L7, Canada 
^Institute for Quantum Computing, University of Waterloo, Waterloo, Ontario, 
N2L 3G1, Canada 

Q_[ ^Department of Chemistry, University of Waterloo, Waterloo, Ontario, N2L 3G1, 

C/3 . Canada 

(Dated: 20 September 2012) 

InAs nanowires with ALJni_ x P or ALdni-^As shells were grown on GaAs substrates 
by the Au-assisted vapour-liquid-solid (VLS) method in a gas source molecular beam 
epitaxy (GS-MBE) system. Core diameters and shell thicknesses were measured 
by transmission electron microscopy (TEM). These measurements were then related 
to selected area diffraction (SAD) patterns to verify either interface coherency or 
relaxation through misfit dislocations. A theoretical strain model is presented to 
determine the critical shell thickness for given core diameters. Zincblende stiffness 
parameters are transformed to their wurtzite counterparts via a well known tensor 
transformation. An energy criterion is then given to determine the shell thickness at 
which coherency is lost and dislocations become favourable. 
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I. INTRODUCTION 



III-V semiconducting nanowires are subject to great interest not only as building blocks 
for solar cells 1 , detectors 2 , and quantum computing devices 3 , but also as novel probes of 
exotic condensed matter phenomena. 4 In all of these applications, optimal device perfor- 
mance depends on the removal of surface states. The elimination of surface states reduces 
unwanted ionized impurity scattering, scattering from rough oxidized surfaces, electron-hole 
recombination, and carrier depletion due to surface traps. There are two methods generally 
used to passivate nanowire surfaces: chemical and structural. Chemical passivation of III- 
V nanowires is typically accomplished with ammonium polysulfide (iVif^S^ 5 ' 6 or organic 
sulfide octadecylthiol (ODT). 7 The effectiveness of chemical passivation declines with time 
so while chemical passivation is a good method to improve ohmic contact formation, it is 
not the best choice for long term passivation of the entire nanowire surface. The alternative 
to chemical passivation is structural passivation where a material with a higher band gap 
is grown around the nanowire core. Improved performance of InAs-InP core-shell nanowire 
field effect transistors (NWFET) has been shown over unpassivated InAs NWFETs. 8 One 
disadvantage of structural passivation is that the shell must be etched selectively to form 
contacts to the core. Fortunately, selective etching is a common practice in III-V device 
processing so etching the shell is generally possible. 9 

The main challenge in realizing core-shell structures is the strain (and dislocations) that 
result from the lattice mismatch. It has been shown that field effect mobility is reduced 
in core-shell structures that have undergone strain relaxation through the formation of 
dislocations. 10 To reduce the likelihood of strain relaxation, a suitable shell material must 
be chosen. In some cases the core and the shell can be lattice matched as in the case of 
GaAs-Alo.52Ino.48P core-shell nanowires, eliminating the possibility of strain relaxation. 11 In 
the case of InAs cores however, there are no III-V materials available (excluding ternary 
antimonides) for lattice matching. As a result, care must be taken not to exceed the thick- 
ness beyond which the formation of dislocations becomes favourable. The traditional way of 
computing the critical thickness of two dimensional epilayers is by the Matthews model. 12 
In recent years there have been a number of models that have addressed the issue of dis- 
locations in the cylindrical core-shell geometry. These models typically have considered 
an isotropic material and calculated the energy required to form a dislocation. 13-19 Some 

2 



models have included the interactions between dislocations and strain fields. 17 More specific 
critical thickness models have been developed for a variety of systems including Si-Ge 20 and 
Al(In)GaN. 21 We have chosen to focus on the In containing ternary alloys AllnP and AlInAs 
for the shell material. InP shells on InAs cores proved to be difficult to achieve via GS-MBE 
so the Al-based ternary alloys were proposed as an alternative. 22 The addition of Al natu- 
rally promotes the formation of a shell by changing the growth mode from Au assisted VLS 
to step flow growth on the nanowire sidewalk Nanowires were grown in a GS-MBE system 
(SVT Associates). Group III species (In, Al) were supplied as monomers from effusion cells 
while group V species (As, P) were supplied as dimers cracked from hydrides. Details of the 
experimental parameters and nanowire characterization are described elsewhere. 22 

II. MODEL 

The strain model presented here was developed using the method of Trammell et al. 
where Si-Ge core-shell nanowires were studied. 20 This model was extended to describe the 
wurtzite structures observed in our InAs-AlInAs(P) nanowires as verified by TEM. In bulk 
form, all of the alloys considered here adopt the cubic zincblende crystal structure. As a 
result, we must convert the relevant parameters for the model from the tabulated zincblende 
values to their wurtzite counterparts. The in-plane equilibrium lattice constant is given by 
a o = Q-zb/Vz and the c-axis equilibrium lattice constant is given by a z Q = ^/8/3 a r Q where 
clzb is the bulk zincblende lattice constant. 23 The well known Martin transformation was 
used to obtain the stiffness tensor with results shown in Table I. 24 

The Martin transformation is motivated by the strong similarity in tetrahedral coordi- 
nation between the zincblende and wurtzite crystal structures. Locally the tetrahedra are 
rotated with respect to each other but otherwise identical. This allows the combination 
of a rotation and an internal strain compensation to estimate the wurtzite stiffness con- 
stants given the zincblende values. The validity of the transformation has been verified by 
experiment in semiconductors that readily crystallize in both the zincblende and wurtzite 
structures such as ZnS. 26 In our study the nanowire cores were InAs, while the shell mate- 
rial was a ternary alloy of ALIni.^As or Al x Ini_ x P. The lattice parameters and the stiffness 
constants for these materials were calculated using Vegard's law. 25 

The nanowire geometry for this model consists of two coaxial cylinders (referred hereafter 
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TABLE I. Non-zero components of the stiffness tensor cy. Shown are the tabulated zincblende 
and calculated wurtzite values. Components of the stiffness tensor for ternary alloys Al^Ini-^P 
and ALjIni-ajAs were computed using Vegard's law. 25 
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InP 


101.1 


56.1 


45.6 


104.95 
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A1P 


188.3 


67.1 


36.9 


164.19 


156.70 


75.41 
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52.22 


44.39 


InAs 


83.29 


45.26 


39.59 


87.83 


110.72 


54.43 


31.54 


13.20 


16.70 


AlAs 


119.9 


57.5 


56.6 


128.87 


153.77 


65.46 


40.57 


26.13 


31.71 




FIG. 1. Geometry of the core-shell nanowire of length L, radius R, core radius R c , and shell 
thickness t in the (r, 0, z) coordinate system. Also shown is the location of a dislocation with 
Burgers vector b. 



4 



as core and shell) of different lattice constant and stiffness tensor. The strain was calculated 
through a two-step equilibration process. In Step 1, the core and shell were considered in 
their initial un-strained state and to be incoherent with each other. The reference state 
was then defined by establishing the condition of epitaxy where a single lattice constant is 
maintained throughout the system. For this model, the reference in-plane and z-axis lattice 
constants were chosen to be those of the core, a ref = a c r and a ref = a c z . In Step 2 the system 
was allowed to elastically relax to the final state. This results in an expression for the total 
strain for k = r,9, z: 



4 = 4 (4 " ml) (1) 
a k 

where m\ = (a k — a c k )/a c k is the misfit strain from the reference state in Step 1, e\ = 
(a j — a c k )l a k i s the displacement strain resulting from Step 2 where a l j is the final lattice 
constant, and i refers to the core or shell. 

The displacement strains from Eq. (1) can be written in terms of the actual displacement 
u r> g jZ in cylindrical coordinates: 27 



du r 

dr 

u r 

£e = — 

r 

du z 



dz 



(2) 
(3) 
(4) 



These displacements were solved by applying a series of boundary conditions appropri- 
ate to the core-shell nanowire geometry. The first boundary condition is a statement of 
mechanical equilibrium due to the absence of external loading: 27 



dz 

where the cr^ are given by Hooke's law (See Reference 28 ): 



1^-^ = (5) 
dUz (6) 
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a r = cne r + c 12 e e + c 12 e z (7) 
ce = c u e r + c u e e + c 12 e z (8) 
<? z = c 12 e r + ci 2 e e + c u e z (9) 

where are the stiffness tensor elements from Table I. Substituting Eq. (7)-(9) into Eq. (5)- 
(6) we obtain the system of four displacement equations; u r and u z in the core and shell: 



d 2 u r 1 du r u 
d 2 u 



dz 2 

Solutions to the displacement equations for the core and shell are: 



3 = o (io) 



(11) 



u r = ar H — (12) 
r 

u z = jz + (j) (13) 

To obtain the eight coefficients a,/?, 7, <j) for both core and shell in the above solution 
we apply boundary conditions arising from geometrical and physical constraints. In the 
nanowire core, as r — > 0, u c r — >■ 00 so /3 C = 0. Continuity at the interfaces was maintained by 
requiring that the lattice constants were continuous across the interface. In the z-direction, 
this eliminates <fi and requires 7 C = 7 s . The interface continuity is further specified by 
maintaining continuity in displacement and stress across the interface: 



u c r (r = R c ) - u s r (r = R c ) = (14) 
a c r (r = R c ) - a s r (r = R c ) = (15) 

Since the nanowire is in mechanical equilibrium, the net force on the surfaces must be 
zero which is expressed as: 

[■R c i-2-k i-B, i-2-k 

\ / a c z rdrd6- \ / a s z rdrd9 = (16) 
Jo Jo Jr c Jo 

a s r ( r = R C ) = (17) 
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The resulting expressions for displacement are obtained through a lengthy solution to the 
system of four algebraic equations, Eq. (14)-(17), to determine the four remaining coefficients 
a c , 7 c ,a s ,/3 s (not shown here due to length). 

The limits of coherency were obtained by comparing the strain energy and the energy 
required to form a dislocation. To determine the strain energy, the displacement solutions 
were substituted back into the expressions for stress in Eq.(7)-(9) and in the expressions 
for strain in Eq.(2)-(4). These expressions were combined to compute the strain energy as 
follows: 



U c = - I dz I dr j dO a c k e c k r (18) 
2 



Jo Jo 

L r R p 2tv 



U s = - I dz I dr I d9 a s k e s k r (19) 



JR C 

where summation over k = r, 9, z is implied. In this model we only consider an edge 
dislocation where an extra plane of atoms is inserted in the shell in the [0001] direction. 
This is the only dislocation we observed experimentally. 22 The formation energy for such a 
dislocation is: 29 

W = Kb 2 a fin ^ - 2J (20) 

where b is the burgers vector of the dislocation, which in this case is equal to the z-axis 
lattice constant of the shell, a is the radius of the dislocation set equal to the core radius 
R c , and r is the cutoff radius that eliminates a mathematical singularity at the core of the 
dislocation. Eq. (20) takes into account the core energy of the dislocation. The cutoff radius 
is set to 6/4 for semiconductors. 30 K is the energy factor and is determined from the elastic 
constants: 



K = (ci 3 + c 13 ) 



C44 (C13 - C13 J 



1/2 

(21) 



.C33 (C13 + C13 + 2c 44 ) 

where C13 = a/ (cnc 33 ). To determine the point at which a dislocation is formed, we deter- 
mine the geometry (core radius R c and shell thickness t) at which the inequality U c +U s > W 
occurs. 
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FIG. 2. Plot of the critical ALjJni-^As shell thickness as a function of InAs core diameter for 
different alloy parameters, x. Also shown are data points from nanowires grown previously denoted 
A, B, and C. 

III. RESULTS 

Figure 2 shows the calculated critical thickness for InAs-Al x Ini_ x As core-shell nanowires. 
As the core radius increases to infinity, the critical thickness becomes constant, meaning the 
model reduces to the thin film case as expected. The results in Figure 2 were compared with 
experimental data. The experimental data points in Figure 3 indicate the shell thickness and 
core diameter measured by HRTEM for a number of InAs-Al x Ini_ x As core-shell nanowires 
with nominal composition x = 0.53 (sample A), 0.36 (sample B) and 0.20 (sample C) as 
reported previously. 22 The measured shell thickness and core diameter was in the range 
of 8-14 nm and 25-40 nm, respectively. In Figure 3, SAD images are shown from InAs- 
Al x Ini_ x As core-shell nanowires. 22 Moving from a)-c) we show a decrease in alloy fraction and 
a subsequent decrease in spot splitting. Spot splitting and spot broadening are associated 
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FIG. 3. SAD pattern of InAs-AL c Ini__ ;c As core-shell nanowires along the [2110] zone axis for a) 
x = 0.53, 3.6% mismatch, b) x = 0.36, 2.5% mismatch, and c) x = 0.20, 1.3% mismatch. A, B, C 
from Figure 2 correspond to a), b) and c) in this figure. 

with strain relaxation, which in this case is due to the epitaxial mismatch between the InAs 
core and the ALJni_ x As shell. Similarly, in Figure 4 we show SAD patterns of nanowires with 
Al x In!„ x .P shells on InAs cores where larger relaxation is evident by further spot splitting. 
For all combinations of shell thickness and core diameter, nanowires in Figure 4 well exceeded 
the calculated critical thickness and will not be discussed further. Instead, we focus on the 
InAs-Al x Jni_ a As core-shell nanowires. 

The shell thickness (8-14 nm) for samples A and B well exceeded the corresponding 
theoretical critical thickness curves for x = 0.53 and x = 0.36 nm. The critical thickness 
in these cases is below 1 nm in the diameter range of 25-40 nm. Hence, the shells in these 
samples are expected to relax as observed in Figure 3. Dislocations were also observed 
directly in HRTEM as described previously. 22 

Experimentally, no dislocations were observed in nanowires with x = 0.20 (sample C) 
consistent with the sharp spots in Figure 3c). We concede however that finding a single 
dislocation in a nanowire that is 1 /xm in length is not an easy task by HRTEM. Energy 
dispersive x-ray analysis (EDX) on these nanowires found that the actual composition of the 
shell in sample C nanowires varied from x = 0.10 to 0.15, somewhat below the nominal value 
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FIG. 4. SAD pattern of InAs-Al^Ini-^P core-shell nanowires along the [2110] zone axis for a) x = 
0.53, 7.1% mismatch, b) x = 0.36, 5.8% mismatch 

of x = 0.20. Due to the proximity of the experimental data for sample C to the x = 0.10 
theoretical curve in Figure 2, we might expect the shells of these samples to be coherent as 
observed in Figure 3. Overall, these data suggest that Al^Ini-^As shells with x < 0.15 and 
thickness below 15 nm on InAs cores (25-40 nm diameter) are coherent and free of misfit 
dislocations. 



IV. SUMMARY 

We have outlined a model to calculate the critical thickness for core-shell InAs-AlInAs(P) 
nanowires. Comparing the strain energy and the dislocation formation energy, we were able 
to determine the point at which the system underwent strain relaxation via the insertion of a 
misfit dislocation. Extending the model by Trammell 20 , and using the transform by Martin 24 
we have taken into account the wurtzite crystal structure of our nanowires. Comparisons 
with SAD patterns of InAs-AlInAs(P) core-shell nanowires with varying Al alloy parameter 
showed different degrees of strain accommodation due to lattice mismatch. The model 
predicted strain relaxation in all samples that showed SAD spot splitting. The Al .2oI n o.8oAs 
shell showed no SAD spot splitting, consisten with the model predictions when the actual 
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Al alloy fraction measured by EDX is considered. 
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